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Thanks to the latest development in the field of Monte Carlo event generators and satel- 
lite programs allowing for a straightforward implementation of any beyond the Standard 
Model theory in those tools, studying the property of any softly-broken supersymmet- 
ric theory is become an easy task. We illustrate this statement in the context of two 
non-minimal supersymmetric theories, namely the Minimal Supersymmetric Standard 
Model with i?-parity violation and the Minimal _R-symmetric Supersymmetric Standard 
Model and choose to probe interaction vertices involving a non-standard color structure 
and the sector of the top quark. We show how to efficiently implement these theories in 
the MATHEMATIOA package FeynRules and use its interfaces to Monte Carlo tools for 
phenomcnological studies. For the latter, we employ the latest version of the MadGraph 
program. 
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1. Introduction 

Many extensions to the Standard Model ol particle physics have been proposed 
over the last forty years. In particular, weak scale supersymmetry, and especially 
its minimal version dubbed the Minimal Supersymmetric Standard Model (MSSM) 
i s one f the most theoretically and experimentally studied alternatives of the 
Standard Model. Supersymmetry is a symmetry relating particles with opposite 
statistics. In the MSSM, one single fermionic (bosonic) superpartner is assigned 
to each of the Standard Model bosonic (fermionic) degree of freedom, with the 
exception of the Higgs sector which contains two doublets of scalar Higgs fields, 
together with their fermionic superpartners. In this way, several conceptual problems 
of the Standard Model are addressed, such as the large hierarchy between the Planck 
and the electroweak scale, gauge coupling unification at high energy or the existence 
of a candidate for dark matter. However, since the superpartners of the Standard 
Model particles have not yet been observed, supersymmetry has to be broken at low 
energies. In addition, in order to remain a viable solution to the hierarchy problem, 
this breaking has to be soft, yielding supersymmetric masses around the TeV scale. 
Therefore, the quest for supersymmetric particles is one of the main topics of the 
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current experimental program at the Large Hadron Collider (LHC) at CERN, which 
is currently exploring the electroweak scale. 

The general purpose LHC experiments, ATLAS and CMS, are currently pushing 
the li mits on the masses of the supersymmetric particles to a higher and higher 
range However, most of the analyses only hold in the context of the so-called 
constrained MSSM, where i?-parity is conserved and where the 105 additional free 
parameters of the most general (i?-parity conserving) version of the MSSM are 
reduced to a set of four parameters and a sign, assuming an organizing principle 
based on unification at high energy. In contrast, there exists a vast variety of (non- 
minimal) supersymmetric models which can be valuable to be investigated both 
from the point of view of the theorist and the one of the experimentalist. 

Studying the hadron-collider phenomenology of these non-minimal models re- 
quires the use of Monte Carlo event generators in order to describe both backgrounds 
and new physics signals. These tools efficiently match a proper modeling of the 
strong interactions (parton showering, fragmentation and hadronization) with the 
calculation of matrix elements underlying the considered hard-scattering processes 
where new physics is expected to appear. For this reason, activities in the field of 
Monte Carlo simulations have been rather intense over the last twenty years, result- 
ing in many advances through the release of automated tree-lev el ma trix-ele ment 
generators, such as Alpg en^, Comix ® , Comp Hep/C alcHep SHE! Hela c ^33, 
MadGraph/MadEventEDMMIQI sherpaQSELZ] anc i Whizard HHES] 

Computing predictions for a given new physics theory with the help of the 
above-mentioned tools requires the implementation of the complete set of associ- 
ated Feynman ru les in these progra ms . This task has been allevi ated by packages 
such as LanHep l 2 ° l 21 l, FeynRules l ^2 | 23 | 24 | 25 | or g ARAH | 26 | 27 | which start from 

the Lagrangian of the theory and export the associated Feynman rules to one or 
several event generators in an automated fashion. However, whereas all these tools 
can address the implementation of the Lagrangian in the usual space-time, only the 
FeynRules packa ge is s uitable to perform ab initio computations within the su- 
perspace formalism l 2 8 | 29 | a na tural framework for supersymmetric model building. 

In a softly broken supersymmetric theory, the field content is embedded into 
chiral and vector supermultiplcts for the matter and gauge sector, respectively, and 
the Lagrangian is written as a sum of four terms. The first two consist in kinetic 
terms for the chiral and vector supermultiplcts of the model. Assuming renormal- 
izability, these terms are entirely fixed by supersymmetry and gauge invariance 
and are hence model-independent. The third piece of the Lagrangian describes the 
interactions among the chiral supermultiplets and is driven by the so-called superpo- 
tential, whilst the last piece contains the supersymmetry-breaking terms. These two 
last parts of the Lagrangian are both model-dependent. In the superspace formal- 
ism, the Lagrangian is written in terms of superfields, in contrast to the component 
field formalism where it is expressed in terms of the usual scalar, fermionic and 
vector fields of particle physics. The main advantage is related to simpler and very 
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compact superfield expressions, in contrast to their component field counterparts. 
This renders the implementation of softly broken supersymmetric Lagrangians in 
high-energy physics tools much more straightforward and bypasses many possible 
sources of mistakes. 

The superspace module of FeynRules contains a set of functions generating au- 
tomatically the kinetic terms of the Lagrangian (expressed in terms of superfields) . 
Subsequently, the implementation of any supersymmetric theory in the FeynRules 
package, and therefore in any of the Monte Carlo event generators interfaced to it, 
consists only in the setting of the superfield content, the model parameters and 
the gauge symmetries of the theory, together with providing the model dependent 
parts of the Lagrangian, i.e., the supersymmetry-breaking Lagrangian and the su- 
perpotential (the latter being expressed in terms of superfields) . Once implemented, 
the Lagrangian is automatically expanded in terms of component fields by Feyn- 
Rules, these objects being the ones handled by the Monte Carlo event generators. 
The Feynman rules are then extracted and the vertices are passed to the different 
interfaced tools. 

This procedure allows for phenomenological investigations of a large class of 
supersymmetric models whose full implementation in a Monte Carlo event genera- 
tors was considered too involved so far. Currentl y, dedi cated interfaces exist to the 
CompHep/CalcHep, FeynArts/FormCalc MadGraph/MadEvent, 
Sherpa and Whizard programs. In addition, a Universal FeynRules Output 
(UFO) of the model can be created. In this case, all the model information is 
exported as a set of Python classes and objects representing particles, parameters 
and vertices. The produced Python library contains all the interactions included in 
the Lagrangian, without any restriction on the allowed Lorentz and/or color struc- 
tures, in contrast to the other interfaces which reject vertices not compliant with the 
(hard-coded) structures supported by the corresponding programs. Presently, the 
UFO format is used by the MadGraph (version 5) and the GoSAM^generatorsfl 

In this work, we present the phenomenological study of collider signatures as- 
sociated to interaction vertices with exotic color structures. We start by describing 
how to implement into FeynRules two examples of supersymmetric theories be- 
yond the constrained MSSM in the most automated possible way, to extract the 
UFO libraries and to use it then with the MadGraph Monte Carlo event generator 
for phenomenological investigations of signatures specific to each model. The choice 
of the MadGraph generator is driven by the color structures which we desire to 
probe, since among the above-mentioned tools, only MadGraph 5 is capable of 
handling non-standard color or Lorentz structures, using, among others, t he f lex- 
ibility of the UFO model format, which is handled internally by Aloha ^H, an 
application automatically writing Helas I 35 | 36 | 37 |38| iib rar i es from Feynman rules 
corresponding to the model under consideration. 



a Even if the GoSam program is also using the UFO as its standard model format, this tool is 
dedicated to next-to-leading order computations, which goes beyond the scope of this work. 
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In our first example, we extend the constrained (i?-parity conserving) MSSM by 
allowing superpotential and supersymmetry-breaking Lagrangian terms in which 
i?-parity is broken ESI. We choose to focus on a set of particular vertices included in 
that model which couple three particles lying in the fundamental representation of 
SU(3) C through a totally antisymmetric tensor of three color (anti)triplet indices. 
For our second example, we turn to phenomenological investigations of scalar par- 
ticles lying in the adjoint representation of the SU (3) c gauge group, as predicted 
by hybrid supersymmetric theories containing representations of both the N = 1 
and N = 2 superalgeb ras I 40 | 41 | 42 | 43 | 44 | 45 ]46| or _R-symmetric supersymmet- 
ric models I 4T | 48 | 49 150| j n those models, the SU(3) C vector supermultiplet of the 
MSSM is supplemented by a chiral supermultiplet, forming hence a complete gauge 
supermultiplet of a N = 2 supersymmetry lying in the adjoint representation of the 
SU(3) C gauge group. This latter supermultiplet contains a vector field (the gluon), 
two (two-component) fermions mixing to a Dirac gluino field (in contrast to the 
Majorana counterpart included in the MSSM) and one scalar particle, commonly 
dubbed as the sgluon field. Even if pairs of sgluons interact with one or two gluons 
through the usual covariant derivative of SU (3) c , hybrid and i?-symmetric super- 
symmetric theories also predict a coupling between a single sgluon and a pair of 
gluons through a dimension-five operator generated by a loop of squarks, the color 
part of the interaction vertex being the symmetric tensor of SU(3) ^4 Since the 
QCD sector is the same in both the i?-symmetric and hybrid supersymmetric mod- 
els, we choose to focus, in the following, only on the i?-symmetric theories for the 
sake of the example. 

Even if we decide to focus on non-minimal softly-broken supersymmetric theo- 
ries predicting interaction vertices with peculiar color structures, a broad class of 
extensions of the MSSM contains additional non-renormalizable interactions in the 
Kahler potential and in the superpotential, leading to non-stan dard Lore ntz struc- 
tures. These new interactions appear, e.g. in the Higgs sector [5 2 | 53 | 54 |^ wnen we 
consider possible new physics effects through higher-dimensional operators. These 
models can be directly implemented into the superspace module of FeynRules and 
consequently exported to Monte Carlo tools for further phenomenological investiga- 
tions in the same way as for the two examples considered in this work. In this case, 
Monte Carlo programs containing a hard-coded set of supported Lorentz structures 
are again not suitable if the relevant structures are too exotic. Once again, the 
MadGraph 5 program, using the UFO and Aloha, allows for the implementation 
of any (renormalizable or not) interaction vertex regardless of the complexity of the 
Lorentz structure. 

This paper is organized as follows. In Section [2j we recall the general features 
of softly-broken supersymmetric Lagrangians and describe the two considered non- 
minimal supersymmetric theories, namely the MSSM with i?-parity violation and 
the minimal i?-symmetric supersymmetric model, and their implementation within 
the superspace module of FeynRules starting from the existing MSSM implemen- 
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tation. In Section |3l we present the extraction of the UFO model format and its 
linking to MadGraph 5, whilst Section [4] focuses on the adopted benchmark sce- 
narios and the description of the high energy physics process of interest. Section [5] 
is dedicated to our phenomenological analyses and our conclusions are finally given 
in Section [6] 

2. Softly-broken supersymmetric theories and their 
implementation in FeynRules 

In this Section, we describe the main features of softly-broken supersymmetric the- 
ories and how to implement the corresponding Lagrangians into the superspace 
module of FeynRules. For further information, we refer to the manual of Feyn- 
Rules 1 ^ and to the details of its superspace module^. We then describe the two 
non-minimal supersymmetric models which we propose to investigate, i.e., the Min- 
imal Supersymmetric Model with i?-parity violation and the Minimal i?-symmetric 
Supersymmetric Standard Model, and finally address their implementation within 
FeynRules. 

2.1. A generic supersymmetric theory in the superspace formalism 
2.1.1. Superfields 

Supers ymmet ric theories are more elegantly formulated within the superspace for- 
malism EHES] Superspace is an extension of the ordinary space-time where we sup- 
plement the usual space-time coordinates x M by the variables 9 and 9, transforming 
as two-component Weyl fermions with opposite chirality (forming thus a Majorana 
spinor (9 a , ()&))■ A superfield <E> is any function of the superspace coordinates and 
can generally expanded as a (finite) Taylor series with respect to the variables 9 
and 9, 



The various coefficients of this expansion are the component fields of the superfield 
$ and form a so-called supermultiplet. The fields z, f, g and d are thus complex 
scalar fields, the spinors £, £, u and p are complex Weyl fermions and is a complex 
vector field. In the expression above, we have introduced the four-vector built from 
the Pauli matrices c M . All the conventions, and in particular those related to spinors 
(invariant products, position of the indices, . . .), strictly follow those presented in 
Ref.[25land Ref.[5l 

The unconstrained superfield $ of Eq. fl} provides a reducible representation 
of the supersymmetric algebra, the number of degrees of freedom embedded in the 
expansion of <f> being too large with respect to those included in the supermultiplets 
representing the N = 1 supersymmetric algebra. Consequently, constraints are im- 
posed on the generic superfield $ so that the number of independent component 



z{x) + 9-i[x) + 9-({x) + 0-0 fix) + 9-9g(x) + O^Ov^x) 
+ 9-9 9-uj(x) + 9-9 9-p(x) + 9-9 9-9d(x) . 



(1) 
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fields is drastically reduced. The two types of superfields necessary to construct most 
of the phenomenologically relevant supersymmetric theories are chiral superfields, 
which include one complex scalar and one two-component fermionic component, 
and vector superfields whose component fields are one real massless vector field and 
the associated Majorana fermion. 

Left and right-handed chiral superfields $l and satisfy the constraints 

D^ L (x,9,9)=0 and D a <P R (x, 9, 9) = , (2) 

where the superderivatives D a and read, in our conventions 153 

Under the constraints of Eq. ^ , the number of degrees of freedom included in the 
expansion of the superfields and § R with respect to the variables 9 and 9 match 
those of the matter supermultiplets. The expansions of these superfields are indeed 
given by 

$ L ( yi 9) = <j>{y) + >/20><<l,{y) - 9-9F(y) , 

(4) 

$ R (y\d) = ^) + V29-^)-9-9F(yl) . 

In these two equations, we have introduced the reduced variable y M = x^ — i9a^9. As 
stated above, the expansion of the superfields and $^ includes hence a complex 
scalar field 0, a Weyl fermion ip and an auxiliary additional scalar field F, which is 
necessary to restore equality between the number of bosonic and fermionic degrees 
of freedom for off-shell supermultiplets. 

Chiral superfields can be declared in FeynRules model files in a way similar 
to the one employed for the declaration of ordinary fields. This extends the original 
FeynArts conventions to the superfield level. Superfields with the same quan- 
tum numbers are collected into classes and each class is implemented as a list of 
Mathematical^ replacement rules. For example, a left-handed superfield 4> could 
be implemented as 

CSF[1] == { 

ClassName -> PHI, 

Chirality -> Left, 

Scalar -> z 

Weyl -> psi, 

Indices -> {Index [Colour] } 

} 

This declares simultaneously a left-handed chiral superfield CSF[1] represented by 
the variable PHI and the associated right-handed chiral superfield $L The latter can 
then be used in FeynRules through the (automatically created) symbol PHIbar. In 



b Mathematica is a registered trademark of Wolfram Research Inc. 
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addition, the implemented superfield carries an index labeled by Colour, illustrating 
the fact that it lies in a non-trivial representation of the QCD gauge group (see 
below). The scalar and fermionic component fields of <J> are represented by the 
symbols z and psi, respectively, and it is important to note that these two fields 
must be independently declared in the FeynRules model file. In contrast, the 
auxiliary _F-field is internally handled and the user does not have to take care of it. 
For instance, the scalar component field could be included in the model file as 

S[l] == { 

ClassName -> z, 

Self Conjugate -> False, 

Indices -> {Index [Colour] } 

> 

We emphasize that the options, such as Indices, which are shared by the component 
field class and the superfield class must be identical. A complete description of the 
particle and superfield classes and properties can be found in Ref. [35] and Ref. 

Vector superfields describe the representations of the N = 1 supersymmet- 
ric algebra containing one massless gauge boson together with the corresponding 
fermionic degree of freedom. These multiplets are embedded into vector superfields 
defined by the reality condition 

V = V 4 . (5) 

However, this last condition leads to numerous non-physical fields which have to 
be integrated out. In a suitable gauge, the Wess-Zumino gauge, these unphysical 
components are eliminated already at the level of the expansion of the superfield V 
in terms of the variables 6 and 0, 



V = Ba^BVf, + i6-0 0-A - iO-6 6-X+ -d-d 6-6D . (6) 



This gauge choice is adopted in the rest of this paper. The component fields included 
in the expansion of the superfield V are a real vector field u M , a Major ana spinor 
(A Q , A") and an additional auxiliary real scalar field D ensuring equal numbers of 
fermionic and bosonic degrees of freedom when going off-shell. 

This superfield can be declared in FeynRules model files in a similar fashion 
as for the chiral superfield. Taking the example of the QCD gauge group, one could 
have 

VSF[1] == { 

ClassName -> GSF, 

GaugeBoson -> G, 

Gaugino -> gow, 

Indices -> {Index [Gluon] } 

> 
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This list of Mathematica replacement rules assigns the symbol GSF to the vector 
superfield VSF [1] and defines the related physical components fields by the variables 
G, representing the gluon, and gow, corresponding to the two-component gluino 
spinor. As for the chiral superfield case, both the gluon and gluino fields must be 
properly declared in the model file and the auxiliary _D-field is treated internally by 
FeynRules. As vector superfields and their component fields are associated to the 
mediation of the gauge interactions, they naturally lie in the adjoint representation 
of the gauge group. This is emphasized by the presence of an index denoted Gluoj^l 
among the properties of the superfield VSF[1], and this index is related to the 
adjoint representation of SU(3) C once this group is defined. This last task could be 
performed through the declaration 

SU3C == 
{ 

Abelian -> False, 

Superfield -> GSF, 

Representations -> { {T, Colour}, {Tb,Colourb} }, 
CouplingConstant -> gs, 
StructureConstant -> f, 
DTerm -> dSUN 

} 

This last Mathematica command declares a variable SU3C representing the SU (3) c 
gauge group and associated to the superfield GSF through the Superfield option 
of the gauge group class. This mapping internally defines the indices carried by the 
superfield GSF as those related to the adjoint representation of the gauge group. In 
addition, two other representations arc declared, the symbols Colour and Colourb 
being the fundamental and anti-fundamental representation indices and T and Tb 
the corresponding representation matrices. The declaration of the gauge group also 
assigns the symbol gs to the strong coupling constant and dSUN and f to the sym- 
metric tensor and antisymmetric structure constants of SU(3), respectively. This 
allows to use consistently build-in FeynRules functions computing automatically, 
e.g., the superfield strength tensor or the chiral superfield kinetic terms. 



2.1.2. Lagrangians 

A generic supersymmetric model is defined by fixing the superfield content and 
the gauge symmetries of the theory. We consider a simple gauge group G and 
denote by g the related gauge coupling constant. We introduce a set of irreducible 
representations 1Z of the corresponding Lie algebra g which are spanned by the 

c Strong interactions play a special role in Monte Carlo programs related to the simulation of 
hadron collisions. Consequently, the names of the quantities related to the QCD gauge group, 
such as the index related to its adjoint representation, follow the strict conventions depicted in 
Ref.|22] 



23,2012 1:10 WSPC/INSTRUCTION FILE paper 



Beyond the Minimal Supersymmetric Standard Model: from theory to phenomenology 9 

hermitian matrices T a , and to each matrix, associate a vector superfield V a . The 
matter sector consists in a set of matter supermultiplets described by left-handed 
chiral superfields $ l lying in the representation 1Z 1 of g, whilst the gauge sector 
is described by a vector superfield V = V a T a . In the case of a semi-simple gauge 
group, i.e., when the gauge group consists in a direct product of several simple gauge 
groups, it is sufficient to generalize the previous setup by associating one vector 
superfield V to each direct factor G J of the gauge group. The chiral superfields of 
the theory then lie in a given representation of the algebras corresponding to each 
of the G J groups. 

In their most general (non-renormalizable) version, softly-broken supersymmet- 
ric Lagrangians can be expressed enti rely with th e help of three fundamental func- 
tions of the superfields of the theory I 56 | 57 | 58 [59|^ namely the Kahler potential K, 
the superpotential W and the gauge kinetic function h a b, given together with the 
soft supersymmetry breaking terms. The Kahler potential is a real function of the 
chiral and anti-chiral superfield content of the theory and contains kinetic terms 
describing its dynamics. The superpotential is an holomorphic function of the chi- 
ral superfields and describes their interactions not driven by the Ncether procedure 
leading to the covariantization of the Lagrangian both with respect to gauge and 
supersymmetric transformations. The gauge kinetic function, symmetric in the two 
adjoint gauge indices a and 6, is also an holomorphic function of the chiral super- 
fields, and addresses the kinetic terms related to the vector supermultiplets of the 
theory. Let us note that the Kahler potential and the gauge kinetic function have 
a very simple form if one assumes renormalizability. In this case, the Lagrangian is 
much simpler and its implementation in Feyn Rules can be automated (see below). 
We however first turn to the most general case. 

Putting everything together, (non-renormalizable) softly-broken supersymmetric 
Lagrangian reads 



where an implicit sum over the chiral and gauge superfields $ and V is understood, 
if relevant, and where W* and h* are the functions conjugate to W and h, respec- 
tively. In this last equation, we have introduced the soft-supersymmetry breaking 
Lagrangian £ so f t and the spinorial superfield strength tensors. In the conventions 




(7) 





a 




where the superderivatives have been denned in Eq. (|3|). 
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The general Lagrangian of Eq. ([7]) can be included in a straightforward manner in 
FeynRules model files using the strengths of the superspace module. One possible 
implementation would proceed in two steps. Firstly, the superpotential, expressed 
in terms of superfields, is stored in the variable SuperPot and the supersymmetry- 
breaking Lagrangian, generally expressed in terms of component fieldfl is imple- 
mented in the variable LSof t. Secondly, for the kinetic part of the Lagrangian, we 
define the quantities 

fi = - K(<S>,&e- 2gV ) + K{e- 29V <S>,&) and p = h ab {<5>)W aa W b a , (9) 

and implement them (in terms of superfields) in the FeynRules model file within 
the variables Omega and rho, respectively. An implementation of the Lagrangian of 
Eq. (J7J is finally given by 

Lag = Theta2Th.etabar2Component [ Omega ] + 
Theta2Component [ rho + SuperPot ] + 
Thetabar2Component [ HC[rho] + HC [SuperPot] ] + 
LSoft; 

In this set of commands, we have introduced the FeynRules functions 
Theta2Component, Thetabar2Component, Theta2Thetabar2Component allowing to 
perform the two- and four-dimensional integrations over the Grassmann variables 
and 0, together with the function HC which transforms an expression into the her- 
mitian conjugate counterpart. Let us note that the left-handed and right-handed 
superfield strength tensors appearing in the definition of p can be cast within the 
functions Superf ieldStrengthL and Superf ieldStrengthR. We refer to Ref. \%5\ 
for more details. 

The Lagrangian density of Eq. {7J still depends on the unphysical F and D 
degrees of freedom of the theory and is expressed in terms of Weyl fermions. In 
contrast, all Monte Carlo event generators which could be used to investigate the 
phenomenology of any softly-broken supersymmetric model require the usage of 
four-component Majorana and Dirac (and not two-component) fermions and physi- 
cal fields. The FeynRules package provides a couple of automated functions allow- 
ing to render the Lagrangian compliant with the requirements of the Monte Carlo 
programs. 

Firstly, the unphysical fields can be integrated out of the Lagrangian by insert- 
ing the solution of their equations of motion. This can be performed by employing 
the functions SolveEqMotionD and SolveEqMotionF for the D- and i* 1 - terms, re- 
spectively, 

Lag = SolveEqMotionD [ Lag ] ; 
Lag = SolveEqMotionF [ Lag ] ; 



For a soft-supcrsymmctry breaking Lagrangian defined in terms of superfields and integration 
over the superspace, we refer to Section 12.31 
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Secondly, the WeylComponents option of the four-component fermion particle 
class allows to link a four-component fermion to its two-component counterparts. Set 
up in this way, FeynRules replaces any occurrence of a product of Weyl spinors by 
the corresponding product of four-component spinors. This task is achieved through 
the automated function WeylToDirac. However, before performing these transfor- 
mations, care must be taken with the gauge group representations in which the 
considered helds lie. Indeed, supersymmetric field theories are built only with the 
help of left-handed chiral superfields. Subsequently, the right-handed parts of the 
four-component fermions are embedded into the corresponding charge-conjugate 
left-handed chiral superfield, lying thus in the complex conjugate representation of 
the gauge group. Mapping Weyl fermions to Dirac and Majorana fermions requires 
then a specific treatment. Taking the example of the QCD gauge group imple- 
mented above, it is sufficient to render the indices Colour and Colourb equal and 
use the properties of the representation matrices T = —T* — — T* to get rid of any 
track of the anti-fundamental representation of SU(3) C , T being the fundamental 
and — T* the anti- fundamental representation matrices, i.e., the complex conjugate 
representation. Hence, a possible implementation would be 

Colourb = Colour; 

Lag = Lag /. { Tb[a_,i_,j_] -> -T[a,j,i] }; 
Lag = WeylToDirac [ Lag ] ; 

For renormalizablc theories, the Kahler potential and the gauge kinetic functions 
take a trivial form 

K ($, $ f ) = and h ab {$>) = S ab , (10) 

and the Lagrangian of Eq. ([7]) simplifies. The complete model is then described by 
the superpotential and the soft-supersymmetry breaking terms alone, the f2 and p 
variables being fully fixed by Eq. pop. In this case, these quantities can be com- 
puted automatically with the help of the CSFKineticTerms and VSFKineticTerms 
functions of FeynRules, the only task left to the user being the implementation of 
the superfield content of the theory, the superpotential and the soft-supersymmetry 
breaking Lagrangian. This renders the implementation of any softly-broken super- 
symmetric theory in FeynRules, and therefore in the Monte Carlo programs inter- 
faced to it, almost entirely automated. The two examples investigated in this work 
follow that approach, the tasks related to programming which are left to the user 
being therefore minimal. 

2.2. The Minimal Supersymmetric Standard Model with R-parity 
violation 

2.2.1. Theoretical framework 

The MSSM is the simplest supersymmet ric m odel, resulting from a direct su- 
persymmetrization of the Standard Model J n the most general form allowed 
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by gauge invariance and renormalizability, the superpotential includes dangerous 
baryon-number (J5) and lepton-number (L) violating interactions predicting, e.g., 
fast proton decays, atomic parity violation or large flavor-changing neutral currents 
1221. To circumvent this issue, a discrete symmetry is imposed, namely i?-parity, 
which forbids all the problematic terms. 

However, neither B nor L conservations are fundamental an( j there is no 
strong theoretical motivation to forbid the i?-parity violating terms of the MSSM 
superpotential. One could rather assume their existence and derive strong experi- 
mental constraints on the corresponding coupling strengths. In general, among all 
the possible i?-parity violating terms, data tends to show a single coupling domi- 
nance effect, i.e., only one single i?-parity coupling can be non- negligible at a time 
I6TI621 

The MSSM is based on the same semi-simple gauge group as the Standard 
Model, Gmssm = SU(3) C x SU(2)l x U(1)y- We associate one vector superfield 
to each of the direct factors of the gauge group. These vector superfields, labeled 
by Vg, Vw and Vy, lie in the adjoint representation of the corresponding simple 
Lie group and are singlets under all the other factors of Gmssm- They are given, 
together with their representations under Gmssm, by 

SU(3) C & Vg = (8,1,0) , 

SU(2) L &V W = (1,3,0) , (11) 

U{1) Y & V b = (1,1,0) . 

The set of physical component fields includes, in addition to the Standard Model 
gauge bosons, the gaugino partners, as shown in Table [1] 

Table 1. The superfields included MSSM gauge sector. The 
gauge boson and gaugino components are given together with 
the representations under the SU (3) c X SU (2)^ xU(1)y gauge 
group. 



Superfield 


Gauge boson 


Gaugino 


Representation 


V B 




B 


(1,1,0) 


Vw 




W 


(1,3,0) 


V G 


9m 


9 


(8,1,0) 



The chiral sector of the theory consists of three generations of Standard Model 
quarks and leptons embedded into three generations of chiral supermultiplets, to- 
gether with their squark and slepton partners, 

Ql = (3,2,i) , W R =(3,l,-h , D i R = (3,l,\) 

~~b ~ ~ 6 ~ ~ 6 ( 12 ) 

4 = (1,2,--) , ^ = (1,1,1) , vj = (1,1,0), 
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where i stands for a generation index and where we have indicated the represen- 
tations of the different superfields under the MSSM gauge group. Let us note that 
since supersymmetric theories must be built only with left-handed (and not right- 
handed) chiral superfields, the associated superfields include the charge-conjugate 
right-handed Standard Model fermions and lie thus in the corresponding complex- 
conjugate representations of Gmssm- For completeness, we have included the right- 
handed neutrino superfields. However, they will be kept non-interacting with any 
other superfield, even at the level of the superpotential. 

In addition, (s)fermion superfields are supplemented with two Higgs chiral super- 
multiplets, embedding thus two SU(2)l doublets of scalar Higgs fields accompanied 
by two doublets of their fermionic Higgsino partners, 

H D = {1,2,~) , Hu = (1,2,±). (13) 

This leads to electroweak symmetry breaking without introducing chiral anomalies 
and to mass generation for both the up-type and down-type fermions. The compo- 
nent fields included in the full matter sector can be found in Table [21 

Table 2. The superfields included in the MSSM chiral sector. The 
scalar and two-component spinor components are given together 
with the representations under the SU(3) C x SU (2)j, X U (l)y gauge 
group. The superscript c denotes charge conjugation. 



Superfield 


Fermion 


Scalar 


Representation 


Qi 






-(S) 


(3,2,|) 




U R 

Jic 
a R 




-it 
U R 

4 


(3,1,-f) 
(3,1,1) 






h 


-(S) 


(1,2,-1) 


Er 


pic 
e R 




e R 


(1,1,1) 


V R 






-it 
V R 


(1,1,0) 


Ho 




H d 


-(3) 


(1,2,-1) 


Hu 








(1,2,1) 
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As stated in Section 12.11 the kinetic and gauge interaction terms of the chiral 
and vector superfields are entirely fixed by gauge invariance and supersymmetry. 
Considering the simplest model with trivial Kahler potential and gauge kinetic 
function (see Eq. (fT0|) ). the kinetic Lagrangian reads, 



k=SU(3) e ,SU(2) L ,U(l) Y 

+ / <\ 2 «<rfi'\> (< ' ■ 2 "" , "- ( 2 '' 1 ' Vl> 

$=chiral content 



(14) 



where the superfield strength tensors have been defined in Eq. ((5]). We have intro- 
duced the non-Abelian vector superfields Vw = V W kT kl Vg = VcT a , where the 
SU(2)l and SU(3) C generators T k and T a are taken, for each term in the sum, 
in the representation appropriate to the considered chiral superfield The three 
gauge coupling constants are denoted by g' , g w and g s . 

Non-gauge interactions among the chiral superfields introduced in Eq. (p~2]) and 
Eq. (|13[) are driven by the superpotential. The superpotential of our model is taken 
as the one of the i?-parity conserving MSSM, denoted by Wrpc, supplemented by 
the so-called A terms 



= Wrpc 



1 ("^) 

2^ijkL l L - L 3 L E R + Kjk^L'Q^L^Rl + 2Kjk tlmn UiuD Rm D Rn . 



In this expression, we have explicitly indicated flavor indices (i, j and k), 
(anti)fundamental color indices (£, m, n) and the SU(2) invariant product (as a 
'•')■ 

Similarly, additional i?-parity violating trilinear scalar interactions derived from 
the form of the superpotential can be included in the soft-supersymmetry breaking 
Lagrangian, 

isoft - L S . RPC l -T ljk V L -V L e k R - T[^l~d k m - \ny mn u Rt d R J k Rn + h.c. , (16) 

where the coupling strengths are included in the T, T' and T" parameters and the 
Lagrangian L s .rpc is the i?-parity conserving MSSM soft-supersymmetry breaking 
Lagrangian. 



2.2.2. Implementation in FeynRules 

In the i?-parity violating scenario described above, the superfield content of the 
theory is identical to the one of the usual i?-parity conserving MSSM. Moreover, 
all particle mixings occurring after electroweak symmetry breaking are also left 
unchanged, since the additional terms do not generate further mixing. Hence, the 



c We do not consider in this work the bilinear Higgs-lepton superfield mixing terms KiL % -Hu. 
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relations linking the gauge- and the mass-eigenstate bases are left unchanged with 
respect to the i?-parity conserving MSSM. Therefore, implementing the considered 
model into FeynRules can be performed very efficiently with a minimal effort. 
Implementing all the modifications in a file labeled rpv.fr, it is then sufficient to 
load this file in FeynRules simultaneously with the built-in MSSM model file, 
denoted by mssm.fr, 

LoadModel["mssm.fr" , "rpv.fr"] ; 

For details on the mssm.fr file, which can be downloaded from the FeynRules 
online model database, we refer to Ref. [25] 

The file rpv . f r includes, on the one hand, the definition of the i?-parity vi- 
olating parameters A, A', A", T, T' and T". This follows the standard rules for 
implementing instances of the particle class (see the FeynRules manual). On the 
other hand, this file contains the implementation of the model Lagrangian. The 
kinetic terms presented in Eq. (fl"4f are directly implemented using the automated 
functions CSFKineticTerms and VSFKineticTerms, 

LagKin = Theta2Thetabar2Component [ CSFKineticTerms [ ] ] + 
Theta2Component [ VSFKineticTerms [ ] ] + 
Thetabar2Component [ VSFKineticTerms [ ] ]; 

The i?-parity conserving MSSM superpotential is stored in the file mssm.fr in 
the variable SPot. The complete Lagrangian containing the interactions driven by 
the superpotential of Eq. (1151) can then be implemented in the rpv.fr file as 

SupW = SPot + 

LLLE[f 1 ,f 2,f 3] Conjugate [PMNS [f 4, fl]] * 

LL[l,f4] LL[2,f2] ER [f 3] + 
LLQD[f4,f5,f3] Conjugate [CKM [f 2 , f 5] ] Conjugate [PMNS [fl ,f 4] ] * 
DR[f3,cl] (LL[l,fl] QL[2,f2,cl] - LL[2,fl] QL[l,f2,cl]) 
1/2 LUDD[f I,f2,f3] Eps[cl,c2,c3] UR[fl, cl] DR[f2,c2] DR[f3,c3]; 
LagW = Theta2Component [ SupW ] + Thetabar2Component [ HC [SupW] ] ; 

In the Mathematica commands above, the symbols LLLE, LLQD and LUDD stand for 
the A, A' and A" parameters of the superpotential, Eps for the fully antisymmetric 
tensor of rank three and LL, ER, QL, UR and DR for the chiral superfields Ll, Er, 
Ql, Ur and Dr. For the first and second terms, the CKM and PMNS matrices are 
explicitly included. This allows to compensate further field redefinitions included in 
the model file (see Ref. [25] for more information) . 

On the same footings, assuming that the soft-supersymmetry breaking La- 
grangian associated to the MSSM with i?-parity conservation is implemented in 
the variable LS, (see mssm.fr), the full supersymmetry-breaking Lagrangian is im- 
plemented from the expression of Eq. (|16|) , 

Tsoft = TLLE[f I,f2,f3] Conjugate [PMNS [f 4, fl] ] * 
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LLs[l,f4] LLs[2,f2] ERs[f3] - 
TLQD[f4,f5,f3] Conjugate [CKM [f 2 , f 5] ] Conjugate [PMNS [fl ,f 4] ] * 

DRs[f3,cl] (LLs[l,fl] QLs[2,f2,cl]-LLs[2,f 1] QLs [1 , f 2 , cl] ) - 
1/2 TUDD[f I,f2,f3] Eps[cl,c2,c3] * 

URs[fl, cl] DRs[f2,c2] DRs[f3,c3]; 
LSoft = LS + Tsoft + HC[Tsoft]; 

where the symbols TLLE, TLQD and TUDD stand for the T, T" and T" parameters 
of the soft-supersymmetry breaking Lagrangian and LLS, ERs, QLs, URs and DRs 
for the scalar component of the Ll, Er, Ql, Ur and Dr superfields, respectively. 
Concerning the presence of the CKM and PMNS matrices, we again refer to Ref. 
[25] 

The complete model Lagrangian is thus given by 

Lag = LagKin + LagW + LSoft; 

In order to render this Lagrangian Lag compliant with the requirements of the 
Monte Carlo programs which are aimed to be used, the auxiliary F- and D-ficlds 
must be integrated out. This can be performed with the help of the SolveEqMotionD 
and SolveEqMotionF commands, 

Lag = SolveEqMotionD [ Lag ] ; 
Lag = SolveEqMotionF [ Lag ] ; 

In addition, all Weyl fermions must be replaced in terms of their four-component 
counterparts, 

Colourb = Colour; 

Lag = Lag /. { Tb [a_ , i_ , j_] ->-T [a, j , i] }; 

Lag = Expandlndices [ Lag , FlavorExpand -> {SU2W, SU2D} ]; 
Lag = WeylToDirac [ Lag ] ; 

The set of commands above first maps the matrices of the anti-fundamental rep- 
resentation of SU(3) C to those related to the fundamental one (see section l2".1.2p . 
and then performs an expansion of all the indices, fundamental or adjoint, related 
to the SU(2)l group. This procedure forces the SU{2)l field rotations from the 
gauge basis to the mass basis, which is necessary for the function WeylToDirac to 
correctly perform the translation to Dirac and Majorana fermions 123. Let us note 
that in the current version of FeynRules, the indices SU2W, SU2D and Colourb are 
non standard. Therefore, the manipulations above have to be performed for any 
model. However, in future versions of the code, the issue of automating this step 
will be addressed. 
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2.3. The Minimal R-symmetric Supersymmetric Standard Model 
2.3.1. Theoretical framework 

The supersymmetry algebra naturally contains a continuous i?-symmetry and there 
is a large class of mecha nisms leading to supersymmetry breaking without breaking 
this i?-symmetry IS3J64] Among the major consequences of the conservation of this 
symmetry, Majorana gaugino masses are forbidden. Therefore, phenomcnologically 
viable models requires the pairing of each gaugino with the fermio nic component of 
a new chiral superfield to form a massive Dirac fermion [ 65 | 66 | 67 | 

The preservation of the i?-symmetry also forbids bilinear Higgs mixing term of 
the superpotential and soft-supersymmetry breaking trilinear scalar interactions. 
Therefore, phenomcnological viability requires the introduction of another set of 
chiral superfields, the i?-partners of the Higgs superfields. Mixing terms between 
the Higgs and these new superfields can be included in the superpotential. Conse- 
quently, this renders the Higgsino fields massive and restores agreement with the 
experimental non-observation of a massless Higgsino field. 

Finally, an unbroken i?-symmetry also ensures that most of the dangerous flavor- 
changing operators which could appear in generic supersymmetric theories are nat- 
urally forbidden. 

Starting from the gauge group Gmssm = SU(3) C x SU(2)l X U(l)y and the 
MSSM superfield content presented in Tables Q] and [2 we add the three chiral 
superfields <J> b , &w and 3>g to the theory. They lie in the adjoint representation of 
the relevant gauge group, 

$ B = (1, 1, 0) , $w = (1, 3, 0) , $ G = (8, 1, 0) . (17) 

Let us note that they form, together with the vector superfields Vb, Vw and Vq, 
three complete vector representations of the N — 2 supersymmetric algebra. 

The Higgsino fields can be rendered massive by mixing the two Higgs chiral 
superfields Hr> and H\j with their i?-partners, the chiral superfields Rd and Rjj, 

i? D = (l,2,i), ifo = (1,2,~). (18) 

The notations related to the component fields of these five new superfields are shown 
in Table El 

Kinetic and gauge interaction terms for the new superfields &w and $g are 
given by the canonical Kahler potential of Eq. (fTU| , 



C K = / dW<9 



$tj$ B + ^ w e- 29 ^ w ^ w + $^e- 29 ^ G $r 



(19) 



where Vw = V W kTk and Vq = Va^T a . The vector superfields V W k and Vg°- have 
been defined in Eq. (fTTj) . and, in this case, the matrices T^ and T a are taken as 
the representation matrices of the SU(2) and SU(3) algebra in the adjoint repre- 
sentation. The Lagrangian of Eq. (Q~9|) contains couplings of the new scalar adjoint 
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Table 3. ij-partners of the vector and Higgs superfields in the 
minimal ij-symmetric supersymmetric model. They are given to- 
gether with their scalar and fermionic components, as well as with 
their representations under the SU(3) C X SU{2)l x U(1)y gauge 
group. 



Superfield 


Fermion 


Scalar 


Representation 


*s 


B' 


&B 


(1,1,0) 




W' 




(1,3,0) 


*G 


g' 




(8,1,0) 


Ru 


*-{%) 


R - ( R °) 
\Ru) 


(1,2,-5) 


Rd 


Mf) 


""{%) 


(1,2, i) 



cr-fields to a single gauge boson and to pairs of gauge bosons through the usual co- 
variant derivatives. The gauge-invariant kinetic terms of the i?-superfields are given 
as in the second term of Eq. (fl4|) . the representation matrices of SU (2) being taken 
in the fundamental representation. 

Interactions among the chiral superfields of the theory is given by the superpo- 
tential 

W(*) = (y u hUkQ 3 L -Hu - (y%^ R Qi-H D iy'' J< n I. II , , 



E 

i=U,D 



(20) 



In the expression above, y u , y d and y 1 denote the 3x3 Yukawa matrices in gen- 
eration space and these terms are also present in the MSSM. In contrast, the other 
terms are related to the i?-partner superfields. The parameters /ij are the (i?-)Higgs 
off-diagonal mass-mixing parameters and the dot products again stand for SU(2) 
invariant products. The A-couplings are the trilinear couplings of the i?-partners of 
the U{1)y and SU(2)l vector superfields with the (i?-)Higgs superfields. 

Soft-supersymmetry breaking originates from hidden sector spurions which pre- 
serve the i?-symmetry. In the most general version of the model, both F-type and 
Z?-type supersymmetry breaking terms are allowed. In the first case, a spurion 
chiral superfield X gets a vacuum expectation value (X) = 9 ■ 9vp whilst in the 
second spinorial field strength tensor W' a gets a vacuum expectation value 

(WD = V29aVD: the quantities vp and vd being related to the order of magnitude 
-^susy of the supersymmetric masses. 

Dirac gaugino masses arise through dimension-five operator generating mixing 
between the fermionic component of the chiral superfield and the gaugino com- 
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ponent of the vector superfield Vu, with k = B, W, G, 



la 



£»ft,i = V ^-ml f d 2 9 W ka $ k . (21) 

This involves a D-type spurion, W k being the superfields strength tensor associated 
to the vector superfield 14 and m\ the corresponding mass parameter. The overall 
normalization factor l/(2g k ) depending on the gauge coupling constant ensures a 
correct normalization of the Dirac mass term. 

The -F-type spurion allows to write down mass terms for the scalar sector of the 
theory, 

£soft,2= £ /dWm|-|^$t$. (22) 

/ iVicTicy 
chiral supcrficlds ouo i 

The sum runs over the chiral superfields included in both Table [5] and Table [3] 
However, both F-type and D-type spurions allow to generate additional mass terms 
for the scalar adjoint tr-fields through the dimension-six operators 

£so ft ,3= £ /dW[(m^-^ + (m3 fe ) 2 -^lTr[$ fe $ fe ]. (23) 



k=B,W,G' 



Assuming real soft-supersymmetry breaking masses, the two real degrees of freedom 
included in each of the complex er-scalar fields do not mix. Consequently, we have 
purely scalar or purely pseudoscalar mass-eigenstates. For the sake of simplicity, 
this hypothesis is adopted in this work and the pseudoscalar state are neglected. 

Finally, the soft-supersymmetry breaking Lagrangian also contains bilinear 
(-R— )Higgs interactions driven by the -F-type spurion, 

xx^ r 

BuHu'Ru + BdHo-Rd + BH u -Ho , 



£soft,4 = / d 2 od 2 e 



M 2 



(24) 



where B, Bjj and Bjy are the corresponding mixing strengths. 

More importantly for hadron-collidcr phenomenology, the new scalar fields also 
couple singly to a pair of quarks and to a pair of gluons through loop-diagrams 
involving squarks, gluinos, neutralinos and charginos. Enhanced by the relative 
magnitude (at the weak scale) of the strong coupling, we only keep the dominant 
interactions involving a sgluon field gq. They can be described by the effective 
Lagrangian 

C eS = a a G [(X L ) f f ^ qf P L T a ^ q ' + (\ R y f ,V qf P R T a *(~ 

+ \ g d abc a Ga G aub G^ c + h.c. . 

In this last equation, we have introduced Vt^ as the Dirac field associated to a quark 
of flavor /, G^b as the gluon field strength tensor and d abc as the SU(3) symmetric 
tensor. The strengths of the interactions included in the effective Lagrangian C e g, 
i.e., the parameters Xl, Xr and X g , are taken as free parameters even if in principle, 
they fully depend on the particle spectrum and sfermion mixings. 



(25) 
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2.3.2. Implementation in FeynRules 

As in Section |2~2| the implementation in FeynRules of the model described above 
can be performed in a very efficient way starting from the existing implementation of 
the MSSM. However, unlike the i?-parity violating model, all the novelties cannot be 
included in an additional file only. The declaration of the gluino and neutralino fields 
must indeed be modified, since they are now Dirac instead of Majorana fermions, 
and this requires a modification of the MSSM implementation. Therefore, we start 
from a copy of the file mssm.fr and include all the modifications presented in Section 

The implementation of the five new chiral superfields as well as the one of all 
the model free parameters strictly follow the rules presented in Ref. [22] and Ref.l25l 
Whilst the implementation of the parameters and the one of the i?-Higgs superfields 
is immediate, let us take the example of the $g superfield to briefly comment on the 
implementation of the i?-partners of the vector superfields of the MSSM. A possible 
set of Mathematica commands for the declaration of the superfield $g reads 

CSF[9] == { 

ClassName -> SGL, 

Chirality -> Left, 

Scalar -> sigG, 

Weyl -> gopw, 

Indices -> { Index [Gluon] } 

} 

where we have assigned the symbols sigG and gopw to the sgluon <jq and the gluino 
g' field, respectively. For the implementation of the component fields, we have, 

S[100] == { 

ClassName -> sigG, 

Unphysical -> True, 

Self Conjugate -> False, 

Indices -> { Index [Gluon] }, 

Definitions -> { sigG[a_] -> sigl[aa] /Sqrt [2] ] > 

>, 

S[101] == { 

ClassName -> sigl, 

Self Conjugate -> True, 

Indices -> { Index [Gluon] }, 

Mass -> Msigl, 

Width -> Wsigl 

}, 

W[100]== { 

ClassName -> gopw, 

Unphysical -> True, 
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Chirality -> Left, 
Self Conjugate -> False, 
Indices -> {Index [Gluon] } 

> 

Before moving on, let us comment on the way used to implement the mixing of 
the real degrees of freedom included in the complex sgluon field. Whilst the symbol 
sigG represents the complex scalar field which is the unphysical component field of 
the chiral superfield $g, the symbol sigl is the real scalar (and not pseudoscalar) 
degree of freedom included in og ■ The pseudoscalar state has indeed been decoupled 
and removed from the mixing relation implemented in the option Definitions of 
the class sigG. The g' field has been associated to the symbol gopw. Denoting by 
gow the gluino g two-component fermion embedded in Vg, the two fields g and g' 
can be related to the (Dirac) four-component gluino field through the option 

WeylComponents -> {gow, gpowbar} 

of the gluino instance of the particle class. 

We now turn to the i?-symmetric Lagrangian presented in Section 12.3.11 All 
the gauge-invariant kinetic terms which are given in Eq. (1141) and Eq. (|19[) are 
automatically computed by issuing the commands 

LagKin = Theta2Thetabar2Component [ CSFKineticTerms [ ] ] + 
Theta2Component [ VSFKineticTerms [ ] ] + 
Thetabar2Component [ VSFKineticTerms [ ] ]; 

The superpotential can also be immediately implemented from Eq. (I20[) , translating 
the textbook expression into a Mathematica declaration of a variable SuperW, 

SuperW = ... 

-luB/2 (HU[1] PhiB RU[2] - HU[2] PhiB RU[1]) + 

ldB/2 (HD[1] PhiB RD [2] - HD [2] PhiB RD[1]) + 

luW PhiW[a] (HU[1] Ta[a,2,i] RU[i] - HU[2] Ta[a,l,i] RU[i]) + 

ldW PhiW[a] (HD[1] Ta[a,2,i] RD [i] - HD [2] Ta[a,l,i] RD [i] ) + 

MUu (HU[1] RU[2] - HU[2] RU[1]) + 

MUd (HD[1] RD[2] - HD [2] RD[1]) ]; 

where the dots stand for the trilinear Yukawa interactions identical as in the MSSM 
and then omitted for brevity ^1 In the implementation above, the SU (2) l con- 
tractions have been expanded, Ta are the symbols representing the fundamental 
representation matrices of SU{2)l, HU, HD, RU and RD are the names of the classes 
associated to the (i?-)Higgs superfields, and PhiB and PhiW to the superfields 
and & w- The superpotential A-parameters of the superpotential are denoted by luB, 
ldB, luW and ldW whilst the /x-parameters are given as MUu and MUd. The associated 
Lagrangian can then be computed as 

Lag = Theta2Component [SuperW] + Thetabar2Component [HC [SuperW] ] ; 
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Finally, the soft-supersymmetry breaking Lagrangian split between Eq. (|2ip. Eq. 
([22|) . Eq. (|23|) and Eq. (l24t can also be implemented directly within the super- 
space module of FeynRules. As an example, we take the gluino Dirac mass term 
implemented in the variable mglno which could be included in the model file as 

mglno = MGl/(2 gs) Ueps[be,al] * 

nc [theta[al] , Superf ieldStrengthL [GSF, be, a], PhiG[a]]]; 
LSoft = Theta2Component [mglno] + Thetabar2Component [HC [mglno] ; 

where MG1 is the symbol associated to the product of the soft supersymmetry- 
breaking mass by the vacuum expectation value of the spurion superfield strength 
tensor W', i.e., mfv£>. For the conventions on the rank-two fully antisymmetric 
tensor UEps and on the nc environment, we refer to Ref. [25l 

3. From FeynRules to MadGraph 
3.1. General features 

Among the whole set of existing automated Monte Carlo tools allowing to address 
phenomenological studies of high-energy physics processes, most all of them con- 
tain restrictions on the color and/or Lorentz structures which can appear in the 
interaction vertices of any model. While all the structures included in the Stan- 
dard Model and the MSSM are generally allowed, vertices with non-standard color 
and/or Lorentz structures are most of the time not supported and must be dis- 
carded from the model implementations. As a consequence, this constrains beyond 
the Standard Model phenomenological explorations which could be performed with 
the help of an automated Monte Carlo tool. One is thus forced to compute squared 
matrix element in a non-automated fashion, e.g., by hand, and implem ent the results 
into non-automated tools such as Herwig | 68 | 69 | Qr p Y thia | 70 | 71| _ Q om g beyond 
two-to-two scattering processes is thus a rather tedious task. 

In this work, we present a way to overco me those limitations with the us e of 
the MadGraph 5 matrix element generator ^ and the UFO model format ^21. 
We focus on the non-trivial color structures appearing in the Lagrangians of the 
theories presented in Section [21 namely the rank-three fully antisymmetric tensor 
in color space e^ mn , where £, m and n are (anti) fundamental color indices, and the 
symmetric tensor of SU(3), d abc , where a, b and c are adjoint color indices. 

The Monte Carlo event generator MadGraph 5 allows for the automated gen- 
eration of tree-level matrix elements associated to any scattering process in a very 
efficient way. In particular, MadGraph 5 is thus a suitable tool to simulate final 
state signatures such as those produced in the high energy proton-proton collisions 
occurring at the LHC. The task of the user consists mainly in specifying the process 
of interest in terms of initial and final state particles, the collision basic informa- 
tion (such as the energy of the colliding beams), a set of basic cuts related to the 
analysis aimed to be performed and, of course, the particle physics model under 
consideration, renormalizable or not. 
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The model library of MadGraph 5 is built upon the FeynRules model 
database. Information is passed from FeynRules to MadGraph through the uni- 
versal output format of FeynRules, dubbed the UFO. This format has been de- 
signed to overcome the above-mentioned restrictions on the allowed Lorentz and 
color structures for the interaction vertices at the level of the Monte Carlo tools. 
The key features are flexibility and modularity through the use of Python classes 
and objects to represent particles, parameters and vertices, saving thus the model 
information in an abstract fashion, also independent from any Monte Carlo pro- 
gram. Consequently, it is universal in the sense that it is not tied to any specific 
matrix element generator. Presently, two generators, MadGraph and GoSam, are 
currently using it as their standard model format, and the future version of Her- 
WIG++ will also be based on the UFO format. 

In order to translate a FeynRules model implementation, such as one of those 
described in Section [21 an automated interface has been included within the public 
version of the FeynRules program. It can be called by typing, in a Mathematica 
session, the command 

WriteUF0[ Lag ] ; 

where the variable Lag contains the Lagrangian of the considered model, expressed 
in terms of the usual scalar fields, (four-component) fcrmions, gauge bosons (to- 
gether with the corresponding ghost fields, if relevant) and tensorial fields of parti- 
cle physics. Let us note that all the states are assumed to be rotated to the mass 
eigenbasis. Once issued, the WriteUFO function internally calls the FeynRules 
core function FeynmanRules in order to compute all the vertices associated to the 
Lagrangian Lag. They are then expanded into a color ® spin basis, 

vai ... anttl ...e n(pu _ . . >pn) = £ C ax...a„ Qij L ^ Hpij . . . >pn) > (26) 

where the variables pi denote the particle momenta and Gy the coupling strengths. 
The quantities C^ 1 '"°" n and L> "' "(pi, . . . ,p n ) are tensors in color and spin space, 
respectively. These tensors act as a basis in the color ® spin space and could be used 
by several vertices, reducing by this way possible redundancies in a model imple- 
mentation. The coupling strengths of a specific vertex are therefore the coordinates 
of the vertex in this basis. Following this structure, vertices are implemented in the 
UFO format by using several different Python objects and their attributes to rep- 
resent the vertices themselves, their Lorentz and color structures and the coupling 
strengths. 

Once all the Lagrangian vertices are decomposed as in Eq. ([26]) . the UFO inter- 
face writes a set of Python files stored in a single directory. These files contain all 
the model information, from the definition of the particle content and the parame- 
ters as instances of the generic UFO particle and parameter class, respectively, to 
the implementation of the (factorized) vertices as instances of the vertex class. 

In order to use the produced UFO model with MadGraph, the outputted files 
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must be copied into the models directory of MadGraph so that they can be used 
for event generation as any other built-in model. 

As soon as a process is specified by the user, MadGraph 5 internally calls 
the Aloha packa ge which generates the subroutines, based on the Helas li- 
brary l -35 | 36 | 3T |38j necessar y to compute the helicity amplitudes related to the pro- 
cess under consideration from the UFO. This allows for an efficient evaluation of 
the associated squared matrix element. Indeed, helicity amplitudes include helicity 
wavefunctions corresponding to specific substructures included in a given Feynman 
diagram, which can be further reused within different diagrams. 

Supersymmetric theories contain, in their most general form, more than sev- 
eral thousands of vertices. A large class of vertices are however vanishing in sce- 
narios of phenomenological interests at colliders. As an example, let us take the 
sector of the up-type squarks (the six scalar partners of the three generations of 
left-handed and right-handed up- type quarks). Since after electroweak symmetry 
breaking, all particles with the same spin, color representation and electric charge 
mix, the six up-type squark gauge-eigenstates undergo a 6 x 6 mixing. Consequently, 
the scalar potential contains 0(1000) four scalar interactions among the up-type 
squark mass-eigenstates which violate flavor most of the time. However, the flavor- 
violating entries of the associated mixing matrix are, in most phenomenologically 
viable benchmark scenarios, vanishing, which renders a large part of the 0(1000) 
interaction vertices equal to zero and irrelevant. 

The presence of zero vertices in UFO model files, or more generally into model 
files for any matrix element generator, considerably slows down event generation, 
since they must be loaded in the computer memory on run time and diagrams 
with a vanishing contribution are effectively generated. Therefore, it may be suit- 
able to remove these vertices from the UFO model files to speed up the evalu- 
ation of the matrix elements, making it much more efficient. This task can be 
done at the FEYNRuLES-level with the help of the WriteRestrictionFile and 
LoadRestriction commands. 

Firstly, the numerical values of all the model parameters, i.e., a full super- 
symmetric spectrum, must be loaded in Feyn Rules. Since the instances of the 
parameter class included in FeynRules models are organized following a structure 
inspired by the Supersymmetry Les Houches Accord (SLHA), the numerical initial- 
ization of the model parameters can be achieved by loading a SLHA file directly 
into the Mathematica session, 

ReadLHAFile [Input -> "susy.dat"]; 

where susy.dat is the filename of the spectrum provided by, e.g., one of the dedi- 
cated existing supersymmetric spectrum generators. The detection of the vanishing 
parameters is performed by issuing the following sequence of Mathematica com- 
mands, 

WriteRestrictionFile [ ]; 
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LoadRestriction ["ZeroValues . rst"] ; 

The WriteRestrictionFile function asks FeynRules to scan over the whole set 
of internal and external parameters and write, in a file named ZeroValues . rst, a 
list of Mathematica replacement rules mapping all the vanishing parameters to 
zero. The LoadRestriction function allows to read this file and load the list of 
rules into FeynRules so that it could be used, at a later stage, by either the UFO 
interface or by any other FeynRules interface to a Monte Carlo program. When 
called, the interfaces apply this set of replacement rules to each vertex derived by 
the function FeynmanRules before translating it to be written to the output files. 
The zero contributions are hence immediately simplified and dropped. If after the 
mapping, a full vertex is numerically evaluated to zero, it is ignored and not written 
at all in the generated output files. 

After this optimization, it can be noted that among all the remaining interac- 
tions, hundreds of vertices consist in four scalar interactions which are, for tree- level 
computations, in general phenomenologically less relevant. Therefore, again for the 
sake of efficiency at the Monte Carlo tool level, it is useful not to include those 
vertices. This task can be done through the Exclude4Scalars option of the Feyn- 
Rules interfaces. In the UFO case, one would have to issue the Mathematica 
command 

WriteUFO [ Lag, Exclude4Scalars -> True ]; 

Let us note that one must keep in mind that the UFO model files generated in 
the optimized way illustrated above are not fully general and highly dependent on 
the considered benchmark scenario (defined here in the file susy.dat), even if at 
the FeynRules level, the model implementation is as general as possible. 



4. Benchmark scenarios and processes of interest 

4.1. The R-parity violating Minimal Supersymmetric Standard 
Model 

The ATLAS and CMS experiments are currently setting impressive limits on the 
masses of the supersymmetric part icles and excluding a significant part of the su- 
persymmetric parameter space j n orc jer to reinterpret the experimental results 
in terms of as many different manifestations of supersymmetry as possible, specific 
benchmark scenarios have been recently proposed ^21 Even if a large emphasis is 
put on the constrained version of the MSSM, the i?-parity violating case is also 
addressed. We therefore adopt as a benchmark scenario for our i?-parity violating 
supersymmetric exploration one of the proposed points, along a line in the i?-parity 
violating MSSM parameter space dubbed the 'RPVS-line 11 ^ 1 . 

This line is inspired by the constrained MSSM. Each benchmark lying on it is 
thus defined by four free parameters defined at the grand unification scale and the 
sign of the bilinear Higgs mixing superpotential /x-parameter. To this restricted set 
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of supersymmetric input parameters, one must supplement the value of one of the 
A" coupling, assuming thus the single coupling dominance hypothesis I ^H^ I anc j 
neglecting any other A", together with all the A, A', T, T' and T" parameters. This 
line is thus perfectly suitable to probe the exotic color structures included in the 
superpotential part of the Lagrangian. 

We fix the universal scalar mass Too = 100 GeV, the universal gaugino mass 
m i/2 — 400 GeV and the universal trilinear coupling Aq — GeV. The ratio of the 
vacuum expectation values of the neutral component of the two Higgs doublets is 
taken as tan/3 = 10, whilst the ^(-parameter is chosen positive. 

For the Standard Model sector, we fix the top quark pole mass to mt = 173.2 
Q e y|73] £ ne bottom quark mass to ra^m^) — 4.2 GeV and the Z-boson mass to 
m z = 91.1876 GeV. The Fermi constant has been taken as G F = 1.16637 x 10~ 5 
GeV -2 , and the strong and electromagnetic coupling constants at the Z-pole as 
a s (mz) = 0.1176 and a(m.z) — 1 = 127.934, according to the Particle Data Group 
Review 

There exists a wide range of experimental measurements constraining the mag- 
nitude of the baryon-nu mber violating A" parameters, such as data related to, e.g., 
K—K systems^ESIZS neutron electric dipole moments rare had ronic -B-decays 
| 78 | 79 |^ nucleo n-antinuc leon oscillations, as well as to double I E0 | 81 | 82 | an( ^ s i n gi e nu _ 
cleon decays IS 3 | 84 | 85 | However, the most restrictive bounds on the A" parameters 
are deduced from cosmological and astrophysical data with the observed flux of 
antiproton cosmic rays^Sl. All these limits are nevertheless not applicable on the 
A"3j7£ couplings, related to the top (s)quarks, for supersymmetric scenarios where 
the lightest supersymmetric particle is lighter than the top quark. In this case, the 
X"3jk parameters are left almost unconstrained, which is particularly appealing for 
LHC physics. We choose A312 = 0.2, since with our setup, at it is presented below, 
the lightest supersymmetric particle is lighter than the top quark. 

The input parameters for the chosen benchmark point are summarized in Table 
21 where we remind that the supersymmetric parameters are defined at the grand 
unification scale. 



Table 4. Input parameters associated to the chosen benchmark scenario in the con- 
text of the i?-parity violating MSSM. All other A"-parameters, together with the A, 
A', T, T" and T" parameters are taken equal to zero. 



m t [GeV] 


m b [GeV] 


m z [GeV] 


Gp [GeV~ 2 ] 




a ( m z) _1 


173.2 


4.2 


91.1876 


1.16637 x 10~ 5 


0.1176 


127.934 


m [GeV] 


m 1/2 [GeV] 


A [GeV] 


tan /9 


sign(M) 


A 312 


100 


400 





10 


> 


0.2 



The soft supersymmetry-breaking masses and interaction strengths at the 
electroweak scale ar e obta ined through renormahzation group running using the 
SPheno 3 package I 87 | 88 |^ which solves the renormahzation group equations nu- 
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Fig. 1. Feynman diagrams associated to ij-parity violating monotop production as generated by 
MadGraph 5. In the notations above, the label nl stands for the lightest neutralino, s~ and d~ 
for strange and down antiquarks and t for a top quark. Attached to the internal lines, stl and 
st2 represent the two top squarks, ssL the left-handed strange squark and sdL~ the left-handed 
down antisquark. 

merically to two-loop order. This program then extracts the particle spectrum and 
mixings at the two-loop level for the Higgs sector and at the one-loop level for all 
the other particles. 

In the sector of the electroweak superpartners, the sleptons are fairly light, with 
masses of 0(200 — 300) GeV, and the neutralino and chargino masses range from 
160 GeV for the lightest neutralino, being the lightest supersymmetric particle, to 
550 GeV for the heavier states. Let us emphasize that since the lightest neutralino is 
lighter than the top quark, the conditions to evade the experimental bounds on the 
A" parameters are fulfilled. The choice of A"3i2 = 0.2 is thus justified. In contrast, 
in the strong sector, the squark masses are a larger, even if still rather modest, 
ranging from 650 GeV to 875 GeV, whilst the gluino is heavier than all the squarks 
with a mass of 900 GeV. 

One can observe that in the scenario depicted above, top squarks and antisquarks 
can be singly produced at the LHC with a sensible rate from an initial associated 
pair of down and strange (anti)quarks. Employing the CTEQ6L1 set of parton 
densities and for a center-of-mass energy of 7 TeV, the corresponding total 
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hadronic cross section indeed reach the level of about 2 pb. Whereas top squarks 
dominantly decay to dijets, their observation in this channel is made difficult by 
the huge QCD background. In contrast, the produced (anti)squark can decay to an 
associated pair of a top quark and a lightest neutralino with a branching ratio of 
15%. 

As i?-parity is violated through baryon-number violating operators, one can ex- 
pect a multijet signature after the decay of both the top quark and the neutralino. 
However, in the considered scenario, the only possible three-body decay of the light- 
est neutralino, via a \"sjk interaction, is not kinematically allowed and the only 
remaining option goes through a four-body decay. As a consequence, the lightest 
neutralino has a long lifetime and decays far outside the detector 1^1. Therefore, the 
missing energy signature, typical from i?-parity conserving supersymmetric scenar- 
ios, is still a convenient observable. This will be investigated in Section [5] where we 
carefully study the i?-parity violating production of a single top quark in as soc iated 
with missing energy, this signature being dubbed as a monotop signature ^U. The 
corresponding Feynman diagrams, as generated by MadGraph 5, are presented in 
Figure Q] 

4.2. The Minimal R-symmetric supersymmetric model 

As stated above, most of the results of the ATLAS and CMS experiments are de- 
rived under the hypothesis of the constrained MSSM. Several reinterpretations exist 
in the context of popular models such as, e.g., the MSSM where supersymmetry 
is broken through gauge interactions, the so-called phenomenological MSSM or the 
Next-to-Minimal Supersymmetric Standard Model. However, there are alternative 
realizations of supersymmetry with highly different properties. Therefore, their in- 
vestigation deserves dedicated studies in order to prepare the reinterpretation of 
the data. The minimal i?-symmetric supersymmetric model is one of such exam- 
ples, since it predicts, among others, the existence of a new color-octet scalar field 
(dubbed the sgluon) which is not present in more conventional experimentally cov- 
ered supersymmetric theories. In this Section, we address the phenomenology of 
such as field, from its production at the LHC to its decay and signature within a 
detector. 

In the model presented in Section 12.31 sgluons can be either singly or pair- 
produced in hadron collisions at high energies from initial states of quarks and 
gluons. Since the sgluon field belongs to the Standard Model sector of supersym- 
metric theorie^l, it is expected to decay mainly into light jets and/or top quarks, if 
kinematically allowed. Let us note even if squarks and gluinos could be lighter than 
the sgluon, which implies that additional decay channels might become possible, 
this case is not considered in this work. Subsequently, the superpartners and their 

f The sgluon field carries a positive ij-parity quantum number, contrary to its gluino g' superpartner 
which has a negative ij-parity. 
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interactions are irrelevant for the phenomenological studies to be performed, and 
we therefore conceive a benchmark scenario where they decouple. 

In order to optimize the generation of the MadGraph model files, the masses 
of the superpartners are set to a very high scale such as, e.g., 1000 TeV, and all the 
mixing matrices related to sfermions, neutralinos and charginos are set to zero. This 
makes all the interaction vertices involving two superpartners or more to vanish, 
so that they are effectively removed from the generated UFO model files by the 
optimization procedure described in Section |3l In addition to the Standard Model 
inputs, it is then enough to fix the sgluon related parameters, i.e., its mass m a 
and its couplings to quarks and gluons Al, Xr and X g defined in Eq. (|25j) . Let us 
note that even if the sgluon mass m a depends, strictly speaking, on several soft 
parameters, i.e., mSj , m?f and m$ G , we follow the simplifying approach to regroup 
all their contributions into a single parameter m a . 

To investigate sgluonic signatures at the LHC, we consider a scenario where the 
sgluon is assumed to have a mass of 500 GeV, a rather collider-friendly value leading 
to a largely visible pair-production cross section at 7 TeV. The sgluon mass is thus 
high enough so that the ti decay channel is open (assuming non- vanishing Xl and/or 
Xr couplings). The production rate of a four top signature is hence enhanced with 
respect to the Standard Model expectation of about 0.3 fb, and the observation of 
such a signal might be a hint of the existence of a sgluon field. We impose that 
the sgluon singly couples only to top quarks and to gluons with strengths given 
by (Al) 3 3 = (Afl) 3 3 = 0.3 and A ff = 1.5 x 10 -4 , all other effective couplings being 
assumed vanishing. 

The numerical values of the free parameters defining our simplified scenario 
are summarized in Table [5J together with the relevant Standard Model inputs. In 
this scenario, the sgluon-pair production cross section reach the level of 0.20 pb, 
whilst the sgluon-pair induced four top production cross section is of about 42 fb, 
i.e., more than about 140 times the Standard Model predictions. The correspond- 
ing Feynman diagrams, as generated by MadGraph 5, are presented in Figure [2] 
Whereas the value of the effective sgluon-gluon-gluon coupling of 1.5 x 10~ 4 is rather 
reduced, it contributes significantly to sgluon-pair production and decays. Indeed, 
the first Feynman diagram of the second line of Fig. [2l involving two effective ver- 
tices, contributes to the sgluon-pair total production cross section by about 15%. 
More importantly, the sgluon branching ratio to a gluon pair cannot be neglected, 
reaching more than 50%, due to the phase-space suppression of the top-antitop pair 
channel. 

Table 5. Input parameters associated to the chosen benchmark scenario in the 
context of investigating sgluon production in the minimal /{-symmetric supersym- 
metric model. We remind that the superpartners are decoupled and thus irrelevant 
and that all the other X{l,r} parameters are vanishing. 



m t [GeV] m z [GeV] a s (mz) 


m CT [GeV] (A L ) 3 3 (\ R fz X g 


173.1 91.1876 0.1176 


500 0.3 0.3 0.00015 
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Fig. 2. Feynman diagrams associated to a four top signature issued from the production of a pair 
of sgluons, as generated by MadGraph 5. In the notations above, the labels t and t~ stand for the 
top quark and antiquark, g for the gluon field and sgl for the sgluon. Let us note that u-channel 
diagrams are omitted since they can be deduced from the corresponding t-channel diagrams (lower 
panel) . 

5. Non-minimal supersymmetric phenomenology 

Event simulation is performed for the LHC collider at a center-of-mass energy of 
^fs = 7 TeV and for an integrated luminosity of 4 fb _1 . Concerning both signal and 
background events, hard scatt erin g matrix elements are described with the Monte 
Carlo generator MadGraph Neglecting all quark masses but the top mass, we 
employ the leading order set of the CTEQ6 parton density fit and identify both 
the renormalization and factorization scales as the transverse mass of the produced 
heavy particles. We then match the events generated by MadGraph with parton 
showering and hadronization as provided by the Pythia program. The version six 
of Pythia ^3 is used for background and sgluon signal events, whilst the version 
8 is used for i?-parity violating monotop signal events due to the exotic color 
structure not compliant with the requirements of Pyt hia 6. We finally perform a 
fast detector simulation with the program Delphes^I, using the publicly available 
CMS detector card. Jets are reconstructed using the anti-fc* algorithm with a radius 
parameter of R — 0.5. The two examples of phenomenological analyses presented 
in this Section are performed with the help of the Mad Analysis 5 package ^1 
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5.1. Monotop production in R-parity violating super symmetry 

The main signatures associated with monotop production can be classified according 
to the top quark decay, 

PP -> t + Xi -> bjj + $ T or bl + $ T , (27) 

where j and b denote light and 6-jets, respectively and £ a charged lepton. The miss- 
ing transverse energy $ T is associated to the lightest neutralino Xi escaping the 
detector invisibly due to its long lifetime as well as to the neutrino appearing in 
the case of leptonically decaying top quarks. Since leptonically decaying monotop 
signatures induce d by i ?-parity violating supersymmetry have been widely investi- 
gated in the past I 94 | 95 |^ we therefore focus on monotop events where the top quark 
decays hadronically. 

The only source of irreducible Standard Model background to an hadronic mono- 
top signal consists in the production of an invisibly decaying Z-boson together with 
at least three jets, one of them being tagged as a b jet. In contrast, there are many 
possible sources of (dominant) instrumental background related to detector effects. 
On the one hand, QCD multijet events with misreconstructed jets produce fake 
missing energy and mimic hence the monotop signature. However, asking for the 
reconstruction of a top quark could help to reject most of these QCD events. On 
the other hand, W plus jets, tt and diboson events where the weak bosons decay to 
nonreconstructed leptons, as well as single top events including non-reconstructed 
or misrecontructed jets might also be a po ssible source of background. 

In recent experimental analyses | 96 | 97 |^ y. b een shown that simple selection 
cuts allow to keep a good control on the background. Inspired, in addition, by the 
parton- level results of Ref. 1911 we start by requiring a large amount of missing trans- 
verse energy $ T > 200 GeV. Within the employed simplified detector simulation, 
this selection cut allows to remove all the QCD multijet events from the background 
and to sensibly to reduce the contribution of the tt channel. We then impose a veto 
on the presence in the final state of any charged lepton (electron or muon) with a 
transverse momentum px > 10 GeV and a pseudorapidity |ry| < 2.5. This last selec- 
tion cut does not affect the signal but sensibly reduce the contributions of Z-boson 
plus jets and M^-boson plus jets events. 

In a second stage, we exploit the presence of a hadronically decaying top quark, 
together with the one of its decay products. We hence demand exactly one 6-tagged 
jet with a transverse momentum px > 50 GeV and a pseudorapidity |ry| < 2.5, 
as well as exactly two light jets with a transverse momentum px > 30 GeV and 
a pseudorapidity |?7| < 2.5. We estimate a 6-tagging efficiency depending on the 
transverse momentum of the jet as presented on Figure 3 (left panel) of Ref. [98l 
together with a charm and light jet mistagging rate (depending also on the four- 
momentum of the jet) as on Figure 6 (right panel) of Ref. 99. Consequently, the 
efficiency of correctly tagging a jet with a transverse momentum of 50 GeV as a 
6-jet is of about 70%, whilst the mistagging rate of a charm (light) jet as a 6-jet is of 
about 40 % (2%). In addition, since the two light jets are issued from a VF-boson, we 



23, 2012 1:10 WSPC/INSTRUCTION FILE paper 



32 Benjamin Fuks 




Fig. 3. We demand a monotop configuration of the final state as described in the text (large 
missing energy, no charged leptons, exactly one 6-tagged jet and two light jets whose the invariant- 
mass is compatible with a W-boson). After applying those selection cuts, we present the invariant- 
mass distribution of the three jets rribjj both for the signal (red) and the dominant sources of 
background, i.e., single top (purple), tt (blue) and diboson (gray) events. 



constrain their invariant-mass rrijj to be compatible with the mass of the VF-boson, 
i.e., rrijj € [mw — 15 GeV, mw + 15 GeV ], with mw — 80.31 GeV. 

The distribution of the invariant-mass rribjj of the three jets is illustrated on 
Figure [3] After applying all the cuts described above, the dominant contributions 
to the Standard Model background are tt plus jets, diboson plus jets and single top 
plus jets events. In contrast, the contributions of all the other sources of background, 
such as VF-boson, Z-boson or QCD events, are reduced to a barely visible level 
and thus under a very good control. We further constrain the system of the three 
selected jets by requiring their invariant-mass rribjj to be compatible with the mass 
of the top quark, lying thus in the range [m t — 20 GeV, m t + 20 GeV ]. The final 
number of selected events is shown in Tabled both for the signal and all the sources 
of background, together with the LHC sensitivity to a monotop signal. We define 
the latter as the number of signal events over the total number of (signal plus 
background) events S/y/S + B, which reaches in our scenario the level of 4.95cr. 

Conversely, a 3er deviation from the Standard Model expectation could be al- 
ready observed for any value of the i?-parity violating parameter A" 312 > 0.11, as- 
suming the supersymmetric spectrum to be unchanged. Since the number of signal 
events is expected to be of the same order of magnitude for moderate superpartner 
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Table 6. Number of remaining events Advents f° r both the different sources of background 
and for the signal after the whole set of selection cuts described in the text. Contributions 
of the VF-boson plus jets, the Z-boson plus jets and the QCD multijet channels are zero 
and thus not indicated. The results correspond to an integrated luminosity of 4 fb~ 1 of 
proton-proton collisions at the LHC collider, running with a center-of-mass energy of 7 
TeV. The LHC sensitivity, defined as the number of signal events over the total number 
of (signal plus background) events S/ Vo + B is also indicated. 



Event sample 


Advents after all the selection cuts presented in the text 


Top-antitop pairs plus jets 


8.2 ± 2.3 


Diboson plus jets 


2.7 ± 0.7 


Single top 


0.9 ± 0.3 


Total background 


11.8 ±2.4 


Monotop signal 


33.2 ± 1.0 


Sensitivity 


4.95 



masses (below or around the TeV scale), the set of selection cuts presented above, 
i.e., standard monotop cuts, would be sufficient to render i?-parity violating super- 
symmetric monotop signatures distinguishable for a large region of the constrained 
MSSM parameter space. 

5.2. Multitop production in minimal R-symmetric supersymmetry 

The top-enriched signature described in Section 14.21 leads to final states with a 
large multiplicity of jets and leptons as the decay products of the four top quarks. 
Therefore, the main sources of Standard Model background is expected to be related 
to rare processes such as the production of four top quarks or the production of a 
top-antitop pair in association with one or several gauge bosons. 

In order to reject a good fraction of the Standard Model background events, 
we require the presence of exactly two charged leptons carrying the same electric 
charge^, with a transverse momentum px > 20 GeV and a pseudorapidity |ry| < 2.5. 
In addition, we require isolation criteria and reject leptons in the case they are at 
a relative distance AR < 0.2 of a jet. 

Moreover, each of the leptons is issued from a leptonically decaying top quark. 
Therefore, it comes accompanied with a neutrino, carrying missing transverse en- 
ergy, and we then ask for a missing transverse energy selection cut, keeping an event 
only if $ T > 40 GeV. 

As stated above, our signal events are rather rich in jets. In particular, we expect 
at least four 6-tagged jets (one for each produced top quark) and four additional light 
jets issue from the hadronically decaying top quarks. Consequently, we only select 
events with at least eight jets, each of them having a transverse energy Et > 20 
GeV, and at least three of the jets are required to be 6-tagged. As in Section 15. 1[ 

g In our simplified detector simulation, the charge of a lepton is always correctly identified, which 
is far from being the case in experiments such as ATLAS or CMS. 
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Fig. 4. Applying selection cuts related to a final state with an important hadronic activity (at 
least eight jets, three b-tagged jets) and containing two leptons with the same electric charge and 
a sensible missing energy, we present the distribution of the Hj< variable defined as the scalar sum 
of the transverse momentum of all the selected leptons and jets both for the signal (red) and the 
dominant source of background, i.e., tt (blue) events. 

our 6-tagging efficiency is estimated as on Figure 3 (left panel) of Ref. [98l together 
with a charm and light jet mistagging rate as on Figure 6 (right panel) of Ref. l99l 
We remind that all the efficiencies depend on the transverse momentum of the jet. 

The important hadronic activity in the final state suggests to investigate the Ht 
variable, defined as the scalar sum of the transverse momentum of all the selected 
leptons and jets. The results are presented in Figure |3J After applying the set of 
cuts described above, the dominant contributions to the background consist in tt 
events, as well as, in a smaller (and negligible) extent, in tt plus one or several 
gauge bosons events. This is also presented in Table [7] where we have omitted the 
non-contributing sources of background. 

Our event selection criteria are very restrictive. However, due to the smallness 
of the signal cross section, these cuts are mandatory to ensure a good background 
rejection. In the case of higher luminosity or more favorable benchmark scenarios, 
one could however lower the requirements by demanding less jets with a transverse 
energy Et > 20, and/or a smaller number of 6-tagged jets. The example of asking 
for six (seven) jets including three (two) 6-tagged jets is shown in the upper (lower) 
panel of Figure [5] In this case, the background clearly dominates over the signal. 

Even if the background rejection is pretty efficient with those cuts, a 4 fb _1 
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Fig. 5. Same as Figure [4] but with a weaker requirement of only six (seven) jets, three (two) of 
them being 6-tagged jets in the upper (lower) panel of the figure. Events related to the production 
of a tt pair in association with one or several additional gauge boson are shown in green, whilst 
diboson events are shown in gray. 
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Tabic 7. Number of events Advents f° r both the different sources of background and 
for the signal after the whole set of selection cuts described in the text. Contributions 
of the single and diboson plus jets, the associated production of a tt pair and one or 
several gauge bosons, the single top plus jets and the QCD multijet channels are zero 
and thus not indicated. The results correspond to an integrated luminosity of 4 fb — 1 
of proton-proton collisions at the LHC collider, running with a center-of-mass energy 
of 7 TeV. 



Event sample 


Advents after the selection cuts presented in the text 


Top-antitop pairs plus jets 


0.5 ± 0.3 


Total background 
Sgluon signal 


0.5 ± 0.3 
0.9 ±0.1 



luminosity is not enough to obtain a sensitivity to the new physics signal. Therefore, 
only lighter sgluons would be accessible at 7 TeV. The results are however promising 
for the 2012 LHC run at a center-of-mass energy of 8 TeV and with an expected 
luminosity of 15 fb _1 . 

6. Conclusions 

In this work, we have illustrated how softly broken supersymmetric quantum field 
theories can be efficiently and easily implemented in the FeynRules package, ex- 
ported then to Monte Carlo event generators and finally be ready for phenomcno- 
logical investigations at hadron colliders. 

Thanks to its superspace module, implementing a supersymmetric theory in 
FeynRules is reduced to the very simple task of providing the superfield and field 
content of the model, defining the free parameters appearing in the Lagrangian as 
well as providing the superpotential, in terms of superfields, and the soft supersym- 
metry breaking Lagrangian. The FeynRules program allows then to automatically 
expand the superfield Lagrangian in terms of the physical component fields, render- 
ing it ready to be exported either to Monte Carlo event generators or to a Python 
library dubbed the UFO model format, containing all the model information. 

The UFO model format is agnostic of any assumption on the new physics model 
under consideration. Any interaction vertex, whatever is the number of incom- 
ing/outgoing particles, the included color and Lorentz structures, can be effectively 
exported, in contrast to the other model format used by automated Monte Carlo 
tools. The latter indeed rely on text files to be parsed and a set of hard-coded 
compliant color and Lorentz structures. 

In this work, we have presented two phenomenological studies of collider signa- 
tures associated to interaction vertices with exotic color structures. Consequently, 
only the MadGraph 5 Monte Carlo generator, using the UFO format, is capable of 
generating events associated to such signatures. We therefore employ it to investi- 
gate monotop production in the framework of the minimal supersymmetric standard 
model with i?-parity violation as well as multitop production in the context of the 
minimal i?-symmetric supersymmetric model. We show that for the chosen bench- 
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mark, the LHC is the perfect machine to unveil the presence of the corresponding 
new physics. 
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